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1.  INTRODUCTION 

Until  recently  the  traditional  theory  of  production  functions  has 
been  steady  state  (static) . Such  models  relate  exogenous  input 
vectors  to  net  output  vectors,  except  for  the  static  Leontief 
input-output  model  which  treats  intermediate  products  explicitly. 
See  [ 1 ] , [ 2 ] . The  limitations  of  a steady  state  correspondence 
between  exogenous  inputs  and  net  outputs  become  apparent  when 
considered  for  processes  like  shipbuilding  where  the  greater  part 
of  production  activity  day  by  day  is  devoted  to  producing  inter- 
mediate products.  The  constant  input  rates  over  some  period  of 
time  (steady  state)  cannot  meaningfully  relate  inputs  to  outputs 
because  the  dynamic  time  histories  of  the  inputs  greatly  affect 
the  result.  Prompted  by  such  problems  for  the  modeling  of  the 
construction  of  ships,  a dynamic  abstract  model  of  production 
functions  has  been  developed  and  reported  on  in  [ ] . This 

abstract  dynamic  model  relates  exogenous  input  and  net  output 
histories  by  mappings  in  function  spaces  with  axioms  which  are 
extensions  of  those  used  for  the  steady  state  theory.  Having 
conceptualized  the  abstract  basis  for  dynamic  production  functions, 
it  is  necessary  to  turn  to  a specific  form  of  a dynamic  production 
function  for  shipbuilding  which  may  serve  production  planning  for 
this  technology.  In  so  doing  one  must  immediately  confront  the 
intermediate  product  outputs  of  the  system,  which  basically 
characterize  the  dynamic  aspects  of  production. 
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The  dynamic  production  function  for  shipbuilding  is  defined  on  a 
directed  graph,  the  nodes  of  which  represent  production  activities, 
as  in  an  activity  analysis,  connected  by  arcs  which  specify  inter- 
mediate product  transfers.  Certain  simplifications  of  the  general 
model  in  [3]  can  be  made  for  shipbuilding.  Starting  with  the 
full  model  so  formulated,  further  specializations  are  made  to 
facilitate  the  development  of  a computer  code  to  express  the 
dynamic  correspondence  between  input  and  output  histories. 
Algorithms  for  this  code  are  given  with  some  examples  of  the 
computations  for  an  aggregated  (few  variables)  form  of  the  pro- 
duction graph  for  shipbuilding. 

The  graph  (network)  of  intermediate  product  production  activities 
has  the  appearance  of  a PERT  NETWORK.  In  the  case  of  a single 
output  in  time,  i.e.  one  unit  and  not  a sequence  of  units,  the 
production  requirements  can  be  generated  by  a backward  calcula- 
tion on  the  network.  However,  for  a sequence  of  output  units  in 
time,  the  model  treats  the  composition  of  single  unit  requirements 
in  a dynamic  mode,  carrying  out  the  interactions  over  time  as 
against  mere  superposition  of  the  single  unit  requirements,  so 
that  criticality  of  an  activity  is  determined  dynamically. 

2.  THE  ABSTRACT  DYNAMIC  PRODUCTION  FUNCTION 

Denote  output  histories  on  a time  interval  [0,+<*>)  by  a vector 
u = (u^,u2,  ...,  u^)  of  such  histories  for  m different  outputs, 
with 

u.  : u.  (t)  = output/unit  time  , E J°'+0°) 
i i i=l,2,  ...,m 

Similarly  denote  input  rate  histories  on  the  interval  [0,+°°) 
by  a vector  x = (xlfX2,  ...,  xfi)  of  such  histories  for  n 
exogenous  inputs.  The  vectors  x and  u are  taken  as  points 
in  product  spaces  of  function  spaces  BM+  of  nonnegative  functions, 
bounded  in  the  Sup  Norm  and  Lebesgue  measurable.  Specifically 
the  Direct  and  Inverse  production  correspondences  are  expressed  as 


BM^ 

x e BM^-^lP(x)  e 2 + 


- 

u e BM’“— =-3L  (u)  e 2 , 


where  P(x)  represents  the  subset  of  vectors  of  output  histories 
realizable  from  the  vector  x of  input  histories,  as  an  element 
of  the  set  of  all  subsets  of  BM^  , and  2L(u)  represents  the 
subset  of  vectors  of  input  histories  yielding  at  least  the  vector 
u of  output  histories,  as  an  element  of  the  set  of  all  subsets 
of  BM^  . In  this  fashion  the  mappings  are  production  functions. 

Thus,  the  abstract  model  is  a function  (correspondence)  between 
input  and  output  histories,  taken  to  satisfy  certain  axioms, 
namely : 

IP.l  P (0)  = {0}  . 

P.2  P (x)  is  totally  bounded  (relatively  compact). 

P.3  P (Xx)  D P (x)  for  X e [l,+«)  . 

P.4  If  ||x||  > 0 , either  P(x)  = (0)  or  there  exists 

u e P(Xx)  for  some  X > 0 and  | |u|  | > 0 . In  the 

second  case  there  exists  for  each  scalar  6 e (0,+®) 

a scalar  X.  e (0,+®)  such  that  (0u)  e P(XC  • x)  . 

o 0 n 

For  each  i e (1,2,  ...,  m}  there  exists  y e BM“  , 

| | y j | > 0 , such  that  P (y)  contains  u e BM^  with 
| |ui| | > 0 . 

P.5  The  correspondence  x -*■  P(x)  is  closed. 

P.5+  The  correspondence  x -*■  P(x)  is  upper  semi-continuous. 

P.6  If  u e P (x)  , ( 0 u ) e P (x)  for  6 e [0,1]  . 

The  norm  of  f e Bm“  as  used  above  is  taken  as 


where  ||f.||  = Sup  (f(t)  | t e [0,+®)}  . These  properties  imply 

similar  ones  for  the  correspondence  u -*■  IL(u)  , to  which  is  added 

JL.O  {x  | x e IL  (u)  , y t 3L(u)  for  y <_  x , | |u|  | > 0}  = : 

Eff  3L(u)  is  totally  bounded. 

3.  THE  GENERAL  ACTIVITY  ANALYSIS  MODEL 


Denote  the  internal  activities  of  production  by  a finite  number  of 
activities  A^,A2,  ...,  A^  , to  which  is  added  an  activity  Aq 
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as  a source  for  exogenous  inputs  and  an  activity  A^  + ^ as  a sink 
for  net  outputs.  Let  x e BM+  denote  a vector  of  input  rate  time 
histories  of  exogenous  factors  available  to  the  production  system; 
including  labor  services,  machine  and  facility  services,  energy 
and  fuels,  and  materials  and  equipment  installed  in  production 
but  not  produced  by  the  system.  A distribution  of  x to  the 
production  activities  A^,A2,  . . . , A^  is  denoted  by  time 
histories  xoi'x02'  ...,  Xq^  • As  a physical  constraint 


Each  activity  A^,A2,  ...,  A^  produces  intermediate  and  possibly 

final  products  (as  spares) , the  time  histories  of  which  are 

denoted  by  V , a = 1,2,  ...,k.  Intermediate  and  final  products 

are  taken  in  a common  output  function  space  BM^  for  time 

histories,  i.e.  V e BM™  , as  histories  of  outputs  per  unit  time. 

a + 

The  possible  transfers  of  intermediate  and  final  products  from 
one  activity  to  another  can  be  displayed  by  a directed  network, 
the  nodes  of  which  represent  production  activities  and  the  arcs 
represent  possible  transfers  of  intermediate  products  as  inputs 
exogenous  to  the  receiving  activity.  Such  a graph  is  illustrated 
below.  The  vector  of  time  histories  of  intermediate  and  final 
product  transfers  from  the  activities  A^Aj,  ...,  A^  , to  an 
activity  A^  (a  = 1,2,  ...,  (k+1) ) is  represented  by 


- I v“  , v”  E EM?  , 
i=l  1 1 


where  v?  is  the  vector  of  time  histories  of  intermediate  products 

produced  by  A.  and  transferred  to  A . The  vector  of  time 
l a 

histories  of  outputs  produced  by  A^  is  represented  by 


k+1 

- I 

j-1 


(oi  1 1 2 f • • • § k)  • 
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Let  u denote  the  vector  of  final  output  histories,  u e . 

Then 


u 


vk+1 

a 


With  the  foregoing  notation,  one  may  express  the  overall  corre- 


spondence x e BM+ 
by 


BM™ 

P(x)  e 2 


under  disposability  of  outputs 


P(x)  = ju  | u = l (vk+1)  , 0 

j v 6v  ' ° v 

( ill  v«) e (*°«  ' J.I  ' {“  ’ 1,2 


e [0,1]  , (v  = 1,2,  . . . , m)  , 


k)  , 


k 


l 

o = l 


x 


oa 


th  th 

where  u^  is  the  v component  of  u , (•)  is  the  v component 

of  (•)  , and  (x0a , va ) -*■  Pa(x0a'va)  denotes  the  production 

correspondence  for  Aa  . The  overall  inverse  correspondence 


u ■*  IL(u)  in  the  same  terms  is: 


I 

IL(u)  = x I xA  = X.  I (xQa)  , X e tl,+»)  , (i  = 1,2,  ...,  n) 

( a = l 

( [ vk+1)  = e u , e e [i,+«) 

\a=l  “ /v 

■ (*c  ■ x *s)  * >•.£ v’)  ■ 

i 


( v = 1,2,  . . . , m) 


(a  = 1,2,  ...,  k) 


where  (V  ) -*•  IL  (V  ) is  the  inverse  correspondence  for  A 
a a a c a 


This  expression  of  the  activity  analysis  model  does  not  explicitly 
allow  for  storage  of  the  intermediate  product  outputs  received 
by  an  activity  for  inputs.  Allowance  for  such  storage  may  be 
handled  in  the  following  way  for  x -*•  P(x)  : 


Introduce 

for 

each  activity 

which  inputs 

IP 

a 

(x  ,v  ) 
oa  a 

histories 

of 

P 

(x  , V ) 

A 

a 

oa  a 

3P  UP  (x  , v ) ) , by  a correspondence 

a a Oa  a r 


BM™  A 

]Pa(xoa'Va)  £ 2 ^ Pa  (Pa  UOa  ' Va  > } = 


| A A /*  ) BMJ 

.v  I V(t)  = J V (t  ) dx  , V e I\(x_,vj;  z 2 

' 0 


a oa 


2)  Let  v^  : v“(t)  = J"v?(-r)dT  , t e [0,+»)  , denote  the 


cumulative  transfers  from  node  A.  to  node  A at  the 

1 a 

time  t , a = 1,2,  ...,  (k+1)  . 


3)  Modify  the  second  constraint  of  P(x)  to: 


/(k+1)  A j \ A 

II  < 1 e ^(^(^-Vj)  , («  = 1,2 k) 


4)  Most  exogenous  inputs  are  not  storable,  being  services  of  one 
kind  or  another.  However  for  those  components  which  are 
storable,  say  {v.,v_,  ...»  v } C { 1 , 2 , ...,  n}  , define 

t 

(x  ) : (x  (t) ) = f (x  (t))  dr  , t e [0,+-) 

oat  v . oa  v . J oa  v . 

1 1 0 1 

and  express  the  last  constraint  as 


£ (X  ) < (X) 

a=l  °a  Vi  = 


i e (1,2,  ...,  s} 


l ) , < (x) 

L . oa  v . = 

a = l l 


i e { s+1 , . . . , n } 


where 


(x) 


V . 

1 


(x(t)  ) 


-f 


(X  (T  ) 


v^dx 


In  the  case  of  the  inverse  correspondence  u + IL(u)  , storage 
of  intermediate  products  may  be  handled  in  the  following  way: 

A 

1)  The  inverse  correspondence  of  IP  (x  ,v  ) -*■  P (IP  (x  ,v  ) ) 

c A a Oa  a a a oa  a 

at  the  follow-on  activity  is  a correspondence 

V 3L  (1L  (V  ) ) which  accumulates  over  time  the  inter- 
a a a a 

mediate  product  inputs  v transferred  to  1L  (V  ) , i.e. 

r a a a 


(x^  , v ) e JL  (3L  (V  )) 
oa  a a a a 


(x  ,v  ) e 3L  (V  ) , 
oa  a a a 


and 


i (1  (V  ))=  j(xn  ,v  ) | V e P (IP  (x  ,v  ))[•  . 
a a a ( oa  a 1 a aa  Oa  a ) 

2)  Modify  the  third  constraint  of  2L(u)  to 


f-  ' Ji  **)  e v“)) ' 


(a  = 1,2, 


• • • / 


k)  . 


8 


3)  If  some  exogenous  inputs  may  be  accumulated,  modify 


- / /kv1  j\\ 

3L  ( 3L  | ' vJ  I to  < 

a\  a\j  = l “// 

third  constraint  as: 


to  accumulate  the  same,  and  express  the 


f(x  ) , • 

\ oa  V1 


(S  ) , U ) 

Oa  v Oa  v , , 

s s + 1 


- l h1  j 

IL  IL  > V J 
a\  a\j  = l a 


(xrt  ) , l 0“ 

oa  v . , x 

n x=l 


Thus  by  a bookkeeping  operation  for  the  correspondence  IP 
and  IL^  , storage  of  inputs  may  be  incorporated . 

A capacity  constraint  for  each  production  activity  A^  needs  to 
be  imposed  for  each  exogenous  input  for  a particular  technology: 


,x  , < ; 1 e I0'+■, 


j e {1,2,  ...,  n}  , stating  that  at  each  point  of  time  the  rate 

of  applying  an  exogenous  input  is  bounded,  denying  infinite  time 

substitution  for  inputs,  i.e.  the  total  input  needed  to  produce 

a unit  cannot  be  applied  in  very  short  intervals  of  time  at 

high  application  rates,  without  limitation.  This  constraint 

applies  particularly  for  the  input  (inverse)  correspondence 

■*  TL^(V^)  . In  the  case  of  the  output  correspondence 

(x„  ,v  ) -*•  IP  (x  ,v  ) the  bounds  imposed  by  the  constraint 
Oa  a a Oa  a r 

k 

7 x < x may  or  may  not  supercede  those  of  (x  ) . < (x  ) . 
a“l  oa  = 1 1 * oa  ] = oa  j 

j e {1,2,  ...,  n}  . 


The  foregoing  general  form  of  the  activity  analysis  model  for 

the  production  correspondences  x -*•  IP(x)  , u -*•  IL(u)  exhibits 

great  complexity  of  structure.  In  the  case  of  x -*-  IP  (x)  , the 
k 

constraint  T x < x permits  a k-dimensional  subset  of  the 
, oa  = 
n a = l 

set  BM+  , of  time  history  distributions  of  exogenous  inputs  to 
the  activities  A^  (a  = 1,2,  ...,  k)  . Each  distribution  yields 
a set  P(Xq^,xQ2,  ...,  xqjc)  of  output  histories  and  IP(x)  is 
the  union  of  such  sets  over  all  permissable  distributions  of 


exogenous  inputs.  The  time  histories  of  total  outputs  of  inter- 
mediate products  at  the  nodes  A^  (a  = 1,2,  ...,  k)  may  have 
similar  freedom  of  distribution  if  shared  by  several  activities 
as  inputs.  Accumulation  of  intermediate  product  and  exogenous 
inputs  adds  a further  time  substitution  possibility  of  similar 
order  of  complexity.  Thus  one  can  see  the  possible  complexity 
of  structure  for  the  overall  correspondence  x -*  3P(x)  . 


In  the  case  of  the  overall  input  correspondence 


3L(u) 


the  net  output  histories  are  required  to  satisfy  only 

l (v  ) = u (v  = 1,2,  . . . , m)  in  the  efficient  case,  which 

i \ a / v v 
a=l 

permits  a potential  latitude  of  distribution  of  sources  of  final 
outputs.  Again  intermediate  products  may  have  a complex  freedom 
of  time  distribution. 


It  is  clear  that  a calculation  in  full  possible  complexity  of 
the  general  activity  analysis  model  is  not  possible.  However, 
in  the  case  of  shipbuilding  certain  simplifications  of  this 
structure  can  be  reasonably  made  as  a practical  matter  in  making 
the  model  specific  to  shipbuilding. 

4 . A SHIPBUILDING  ACTIVITY  ANALYSIS  MODEL 


For  this  specific  case  it  is  reasonable  to  assume  that  each 
activity  A^  (a  = 1,2,  ...,  (k-1))  yields  only  a single  inter- 

mediate product,  with  A^  yielding  completed  units.  Each 
activity  can  be  modeled  in  terms  of  time  variable  input-output 
coefficients  driven  by  a scalar  intensity  function  representing 
the  time  history  of  operation.  The  resulting  system  is  an  open 
Leontief-like  dynamic  activity  analysis  model.  Beyond  the 
intrinsic  bounds  (xQc[)j  ( j = 1 , 2 , . . . , n)  on  the  rate  of 

applying  inputs  to  an  activity  A^  , there  may  be  certain  fixed 
facilities  and  major  equipment  (cranes)  which  limit  further 
the  intensity  with  which  the  production  activity  Aq  may  be 
carried  out.  Such  bounds  determine  the  minimum  time  for  completion 
of  a unit  of  output  at  A^  . These  minimal  times  set  lower 
bounds  to  the  time  lags  encountered  in  production.  The  actual 
lags  are  determined  endogenously  by  the  intensity  function  used 
in  operating  the  production  activity. 
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A simplified  example  of  the  expression  of  shipbuilding  in  terms 


of  activities 
the  accompanying  figure. 


A1 ,A2' 


with 


•8  > \+1  = Ag  , is  given  in 

The  construction  activity  for  erection 


of  hull  structure  is  denoted  by  A^  , which  yields  a hull  structure 
as  an  intermediate  product  input  to  activities  A2  , A^  , A^  , A^ 
of  installing  the  propulsion  unit,  the  electric  plant,  the 
command  control  and  auxiliary  systems  respectively.  The  activities 
A2  , A-j  , A ^ yield  installed  propulsion  unit,  installed  electric 
plant  and  installed  command  control  system,  respectively,  as 


intermediate  products  to  activities  Ac  and 


yielding  ship 


equipment  and  outfitting  and  furnishings  respectively  which  deliver 


the  same  as  intermediate  products  to  an  activity 
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of  construc- 


tion servicing  which  yields  the  completed  unit  to  the  sink 


At  the  aggregate  level  in  which  the  example  is  expressed,  A^, 
delivers  the  same  unit  to  A2  , A^  , A^  and  A^  simultaneously, 
and  A2  , A-j  , A^  each  deliver  the  same  completed  unit  simul- 
taneously to  Ag  and  A^  , with  the  latter  delivering  the  same 
completed  unit  to  A„  . For  bookkeeping  purposes  with  our  notation, 

O 

one  may  have  A^  deliver  simultaneously  one  fourth  of  a completed 
unit  to  A2  , A-j  , A ^ and  A^  with  the  latter  delivering  simul- 
taneously one  half  of  a completed  unit  to  Ag  and  A 7 , which 
in  turn  deliver  together  a single  completed  unit  to  Ag  . At  a 
more  disaggregated  level  of  modeling,  the  activities  can  be 
arranged  perhaps  so  that  partitioning  of  output  is  not  required. 


Returning  to  the  general  expression  of  our  model,  let 


Z : Z (t) 
a a 


t c [ 0 , +»)  , a e {1,2,  ...,  k} 


denote  a scalar  intensity  function  for  the  production  activity  A 

Ct 

This  intensity  function  drives  input  coefficient  functions 


II 

0 

s* 

lla  l 

(t) ,a° (t) , 

• • * an (t)  1 1 

, a = 1,2, 

m2  = 

II  sj 

(t ) ,a*(t) , 

••• 

II  , o=l,2 

for  exogenous  inputs  and  intermediate  product  inputs  respectively. 

The  coefficient  a1  (t)  represents  at  the  time  t the  amount  of 
th  ^ 

the  j system  exogenous  input  per  unit  time  required  per  unit 

intensity  of  operating  the  activity  A . The  coefficient  a?(t) 

& 1 


aq  source  of  exogenous  inputs 

A1  HULL  STRUCTURE 
a2  FROPULSION  UNIT 

a3  electric  plant 
a4  command-control 
a5  auxiliary  systems 
A6  ship  equipment 
A7  outfitting  and  furnishings 

Ag  CONSTRUCTION  SERVICING 
Ag  SINK  OF  FINAL  OUTPUTS 


Va  £ BMf  ' Va  = (V1 ' V2 Vk-l)  C BM+' 


FIGURE  1 
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represents  at  the  time  t the  amount  of  the  activity  output 

per  unit  time  required  per  unit  intensity  of  operating  the  activity 


These  coefficients  are  taken  time  variable  to  allow  for  "learning 
effects,"  and  even  possibly  to  allow  for  technological  change. 

They  are  nonnegative  with  no  restriction  otherwise  imposed  on  the 
coefficients  . However  for  it  is  assumed  that: 

(a)  There  exists  for  each  t e [0,+»)  and  a e (1,2,  . ..,  k}  , 
a positive  coefficient  a ^ ( t ) for  at  least  one 

j e (1,2,  ...,  n}  . 

(b)  For  each  system  exogenous  input  j e (1,2,  ...»  n}  , 

there  exists  a time  t e [0,+®)  and  at  least  one  production 
activity  , a e {1,2,  . ..,  k}  such  that  a“ (t) 

is  positive. 

Assumptions  (a)  and  (b)  state  that  an  exogenous  input  is  essential 
for  production  and  all  exogenous  inputs  are  required  somewhere 
at  some  time. 

The  output  coefficients  driven  by  the  intensity  function  are 
taken  as 


ca  = I |ca(t) 


, c (t)  > 0 , a = 1 , 2 , ...,k. 


They  are  understood  to  be  the  fraction  of  a unit  completed  per 
unit  time,  at  the  time  t , per  unit  intensity  of  operating  A^  , 
if  the  output  is  a sizeable  discrete  unit,  otherwise  they 
represent  amount  (in  bulk) . The  same  convention  will  be  applied 
for  the  intermediate  product  inputs.  In  this  way  the  expression 
of  the  model  may  be  made  in  simpler  terms  than  that  required 
for  handling  large  discrete  units. 

The  bounds  (xQa)j  on  the  time  histories  of  the  rates  of  applying 


system  exogenous  inputs  to  the  activity 
bound  Z : Z (t)  , t e [0,+,»)  on  the 


imply  a time  history 


Z^  : Z^  (t)  , t e [ 0 , +«> ) on  the  intensity  function  given  by 


Z (t)  = Min 
a . , 


j c {1,2, 


(a  = 1,2,  . . . , k) 
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otherwise  an  exogenous  input  may  be  applied  at  a higher  rate  than 
is  physically  possible.  In  addition  certain  exogenous  inputs  may 
be  bounded  as  to  rate  of  application,  because  of  fixed  facilities 
and  major  equipment  for  a given  shipyard.  Hence  in  modeling 
the  production  system  of  a given  shipyard  there  may  be  additional 
constraints  on  the  rates  at  which  certain  exogenous  inputs  may 
be  applied,  say 


o = l 


(xoa(t))j  = (xo(t))j 


t e [ 0 , +<*>) 


for  j e { v ^ , v 2 , ...,  vp}  c {1,2,  ...,  n} 


In  the  case  where 


the  exogenous  inputs  v]/v2'  •••»  VF  each  apply  to  single 

distinct  activities  •••»  , ^ui'w2'  C *1,2,  •••» 

such  additional  constraints  may  be  merged  with  the  bounds  (x  (t) ) . 

_ _ Oa  j 

to  obtain  a revised  set  of  bounds  Z : Z (t)  , t e [0,+°°)  , 

a a 

Vj  , j e {1,2,  ...,  F} 


(a  = 1,2,  ...,  k)  . When  an  input 
is  shared  by  two  or  more  activities,  say  (a  = a.,a~,  ...,  a_)  , 


an  allocation 
where 


, (i  — 1,2,  ...,  I)  , I e {1,2,  ...,  R}  , 


l ( X { t)\  = (X  (t)) 

i=iioai  /vj  ° 


as  a distribution  of  resources,  is  only  one  of  many  possibilities. 
The  full  possibilities  in  production  involve  all  the  alternative 
bounds  on  the  intensity  functions  so  generable,  which  permit 
important  factor  and  time  substitution  possibilities.  By  this 
supplementary  imposition  of  bounds  on  the  intensity  functions 
Z , the  exogenous  inputs  v1fv_,  ...,  v are  deleted  from  the 
vector  x and  included  as  part  of  the  physical  structure  of  the 
specific  shipyard  being  modeled. 


The  minimal  delays  or  time  lags  T^  in  production  for  the 
activities  A^,A2,  ...,  A^  are  determined  from  the  bounds  by 


T 

a 


(t)dT  = total  amount  of  input  j 
required  per  unit  output, 


where  ja  is  the  restricting  system  exogenous  input  determining 
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If  the  activity  A is  operated  below  the  bound  Z , i.e.  it  is 

a a 

not  limiting,  the  time  lag  encountered  depends  upon  the  time 
history  Z (t)  , t e [0,+«)  at  which  A is  operated.  Thus  the 

OL  OL 

time  lags  in  production  are  determined  endogenously  as  follows: 


If  t is  the  time  at  which  a unit  enters  production  in  the 
activity  A^  , and 
is  completed  under 
determined  by 


Vfc) 

Z (t) 
a 


is  the  increment  of  time  at  which  it 
, t e (0,+«»)  , the  lag  (t)  is 


t+Vfc) 


/ 


t ( T ) Z (T  )dr 


= total  amount  of  input  j 
required  per  unit  output. 


One  final  simplification:  Since  each  activity  A^,A2»  ...»  A^_^ 

produces  a single  intermediate  product  output,  the  vector  v^ 
will  be  taken  as: 


va  * • • , ^ ^ — 1 r 2 , ...,  k)  • 

Turning  now  to  the  output  correspondence  (x  ,v  ) -*■  P (x  ,v  ) 

r oa  a a oa  a 

for  A , since  fixed  time  lags  are  not  used,  calculations  for 
a 

this  correspondence  may  be  made  discretely  on  a time  grid 
t = 0,1,2,  ...,  N,  ...  for  some  convenient  unit  of  time. 

The  intensity  function  Z^  is  constrained  to  satisfy  (a  =1,2,  ...,  k) 

(1)  Z (t)  > 0 , t = 0,1,2, 

a = 

(2)  Z (t)a“(t)  < (x  (t))  . 

a 3 = oa  j 

t 

(3)  Z (t)i“(t)  1 l v“(t) 

01  1 " T = 1 X 

Z (0)aa(0)  < 0 

a l = 

(4)  Z (t)  < Z (t)  , t = 0,1,2,  N,  ... 

a = a 

Conditions  (1)  and  (4)  merely  reflect  lower  and  upper  bounds  on 
the  intensity  function.  Condition  (2)  reflects  a given  distribution 
of  exogenous  inputs  to  the  production  activities,  and  requires 
that  these  allocations  be  not  exceeded.  The  full  possibilities 
in  production  are  generated  by  the  union  of  the  results  for  all 


«•*,  N,  ... 

(j  = 1,2,  ...,  n 
' (t  = 0,1,2,  . . . , N,  . . 

(t"1*  it-  = 

l Z (t  ) a?  (t)  | _ 

in  a 1 " 


T = 0 


1,2, 

1,2, 


N,  . . 
(k-1) 
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such  distributions.  Constraint  (3)  allows  accumulation  of  trans- 
ferred intermediate  product  inputs  and  requires  that  the  current 
input  in  a unit  of  time  not  exceed  the  stock  available. 

These  constraints  are  expressed  with  the  convention  that  inputs 
are  applied  at  uniform  rates  during  each  unit  interval  and  are 
charged  at  the  beginning  of  each  period  for  the  ensuing  unit  of 
time,  with  output  assigned  to  the  end  of  the  period. 

The  maximal  value  of  (t)  for  t = 0,1,2,  ...,  N,  ...  permitted 
by  these  constraints  may  be  calculated  as  follows: 

From  (2) 


(x  (t)  ) . 

Za(t)  < Ra(t)  , Ra(t)  :=  Min  J — 22 1 

j eS° ( t ) ( 

Sa(t):  =|j  | j e {1,2,  ...,  n}  , a“(t)  > 0 | , 
t — 0,1,2,  ...,  N , ... 


Note : 
a“  (t) 
of  Sa 


Assumption  (a)  for  flA  implies  Sa (t)  not  empty.  If 
varies  in  time  only  by  "learning  effects,"  the  dependence 
upon  t is  moot,  but  not  for  technological  change. 


Let  £a  (t)  : = i | i e (1,2,  ...,  (k-1) } , aa  (t)  > o|  . Define 


Wa(0) 


| Min  ^Ra  ( 0 ) , Za(0)J 


for  £a(0)  = 0 
for  la(0)  / 0 . 


Take  Z (0)  = Wa(0)  , and  from  (3) 
a 

Z (l)aa(l)  < va(l)  - Wa  (0)  • Sa(0)  . 

0 1=1  l 


Za(l)  < Wa(l)  where 


Then  (2),  (3),  (4)  imply 
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Wa(l)  = Min  |V*(1)  , Za(l)  , Na(l)j 


1 1 

l va(T)  - Wa(0)  • i“(0) 

IZO , I“(l)  * 0 

i“(l) 


Nu  (1 ) -< 


ie[a(l) 


+ « , 


r<  i)  - » 


Following  the  policy  of  taking  the  maximal  value  for  Zq  , take 

Z (1)  = Wa(l)  . 
a 

Continuing  in  this  fashion,  the  general  solution  is  found  to  be: 


where 


Z (t)  = W (t)  , t = 0,1,2,  ... 
a 


(Min  [R°(0,  » 2 (0)]  for  £a(0)  = 0 
Wa(0)  - L a J 

( 0 for  l (0)  0 

Wa(t)  : = Min  [r°  (t)  , Zq  (N)  , Na(t)j  , t = 1,2,3,  ... 

I v“(t)  - l Wa  (x  ) a?  (t  ) 

Min  1=2 , Ia<t)  * 0 

|i«J-(t)  ifli(t)  J 


N (t)  : = 


+ 00  , 


P< t)  = 0 


t = 1,2,3,  ... 


Thus  given  a vector  v^  = (v^v^,  ...,  vk-i)  of  time  histories 
of  intermediate  product  transfers  to  the  activity  A^  , the  time 
history  of  output  from  A^  is  given  by 


V (t)  = c°(t)  • W°(t)  t = 0,1,2,  ...  . 

a 


In  the  case  of  a shipyard,  the  transfers  v?  (i  = 1,2,  ...»  (k— 1 ) ) 
are  likely  to  be  positive  for  only  a few  indices  i , and,  in 
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most  cases,  positive  for  only  one  transfer.  See  the  example  at 
the  beginning  of  this  section.  In  substance  a single  time  history 
of  transfer  is  made  from  each  node.  Nodes  A2  , , A^  , A^ 

receive  a single  input,  nodes  Ag  and  A^  receive  three  inputs 
and  node  Ag  receives  three.  However  this  may  be,  when  more 
than  one  transfer  is  made  from  a node,  not  in  fixed  proportions, 
the  total  output  possibilities  for  the  entire  system  encompass 
all  such  possible  distributions. 


The  foregoing  solution  has  been  based  on  nonaccumulation  for  all 
exogenous  inputs.  If  inputs  vi'v2'  va  can  accumulated, 

then  for  j e {v.,v_,  ...,  v } constraint  (2)  should  be  expressed  as 

i ^ cl 

t 

(2a)  l Z (t ) aa  (t)  < (S  (t) ) ; i = {1,2,  ...,  a} 

T = 0 ° Vi  = °“  Vi 


for  t = 0,1,2,  ...  . With  this  modification  the  general  solution 
becomes : 

Za(t)  = Wa(t)  , t = 0,1,2,  ... 

where 

(Min  ["r0  ( 0 ) , 2 (0)  , Ra  (0)1  for  7“  (0)  = 0 
W (0)  = J L a -I 

( 0 for  £a(0)  jt  0 

Wa(t)  : = Min  |V(t)  , Ra(t)  , Zq  (t)  , Na  (t)J  , t = 1,2,  ... 

Sa (t)  :=  {i  I i e {1,2,  ...,  a)  , a®  (t)  > oj 
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Ra (0):  =( 


Min 

ieSa (0) 


(xo  (0)) 

Oa  v ^ 

aa  (0) 
vi 


, Sa(0)  ? 0 


+ ® 


sa(0)  = 0 


N° (t): 


Min 
ie[a (t) 


l v?(t)  - l Wa(x)a“(T) 

_ n 1 n 1 


T — 0 


T=0 


a“  (t) 


, I°(t)  ft  0 


t = 1,2, 


+ oo 


r ( t)  = ^ 


V (t)  = ca  (t) Wa  (t)  , t = 0,1,2,  . 
a 


The  inverse  correspondence  for  a production  activit 

ft 

expressed  in  cumulative  output  required,  say  V 


A 

a 

(t) 


may  be 


9 

Let 


a a 

t c [0,+«>)  . A finite  horizon  T will  be  considered, 
t = T , T-l , T-2 , 

Let  N 

a 

= 0,1,2,  ...  the  constraints  on  the  in- 
are  : 


. , denote  the  number  of  time  units  preceding 
be  the  minimal  o such  that  V(a)  = 0 . In  terms  of 


= T 
T 


the  time  variable 
tensity  function 


Za(o> 


(1) ' 

Z ( ct ) >0;  o=l,2,  ... 
a = 

(2)  * 

(xoa(a,)j  = ; j 

— 1,2,  • • • t n , 

(3)  ' 

0°(o)  > 0?  (1)  - l Z.  (T 

)ia(t)  ; 1 = 

1 = 1 T=1  1 

l o=2,3 

(4)  ’ 

Z (a)  < Z (a)  ; a = 1 , 2 , . 

a = a 

a A ft  , 

• • 

(5)  • 

ca (1)  Z (1)  < V (0)  - V (1) 
a = a a 

„ (o-l) 

c (a)  Z (a)  + T ca  (t)  Z 

A ft 

(t)  < V (0)  - V 

a a 

T = 1 

= a a 

1,2,  ... 

, (k-1)  , 


Constraints  (1)'  and  (4)'  have  the  same  meaning  as  (1)  and  (4) 
for  the  output  correspondence.  The  intensity  function  is  driven 
by  the  constraint  (5)'  taken  with  (4)'.  A policy  of  producing 
as  much  as  possible  for  current  requirements  is  involved  in  using 
(5)',  subject  to  the  capacity  constraint  (4)'.  Constraints  (2)' 
and  (3) ' drive  exogenous  input  requirements  noncumulatively  and 
intermediate  product  transfers  cumulatively. 
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Consider  a = 1 . Then 


Z (1)  < 
a = 


V (0)  - V (1) 


a , , , 
c (1) 


, Z (1)  < Z (1) 
a = a 


Z (1)  < Min 
a = 


V (0)  - V (1) 

— 2 , z (1)  =:Ea(l) 

ca(l)  a 


Take  Z^d)  = E d)  in  order  to  produce  as  much  as  is  possible, 


For  o = 


A 

V (0) 

Z (2)  < -2 

a = 


V (2)  - ca(l)Ea(l) 
a 

ca  (2) 


, Z (2)  < Z (2)  , 
a = a 


V (0) 

Z (2)  < Min  — 


V (2)  - ca  (1)  Ea  (1) 
a 

ca  (2) 


' V2) 


Take  Z (2)  = E (2)  and  continue  in  this  fashion  to  obtain; 
a 


Za(o)  = E (o ) , 0 = 1,2, 


where 


Ea(a)  : 


* (0-1) 

v CO)  - V (o)  - l c“(T)Ea(T) 
= Min  a a T=1 


c°  (0) 


' Vo) 


0 2,3,  ••• 


„ V (0)  - (1) 

E (1)  = Min  — , Z (1) 

ca(l)  0 


1 
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Then  the  noncuraulative  exogenous  inputs  required  are  taken  con- 
servatively as 

(xoa(°))j  = a“(o)  • Ea(a)  , c = 1,2,  ...  , j = 1,2,  ...,  n 

and  the  cumulative  intermediate  product  transfers  required  are 
taken  conservatively  as: 


Aa  , . Ai  , 

vi  (°)  = v.  (1) 


( a — 1 ) 

l i“(T)Ea(T) 

til  1 


0 — 2,3, 

1 = 1,2, 


(k-1) 


where  v^(l)  is  a terminal  value  equal  to  the  total  intermediate 
product  transfers  required  for  Va(0)  . 


The  foregoing  recursive  formulas  for  the  intensity  functions  for 

A A A 

the  correspondences  (x  , v ) =>  ]P  (x  , v ) , V -*•  jL  (V  ) , 

oa  a a oa  a a a a 

a = 1,2,  ...,  k enable  the  calculation  of  the  overall  corre- 
spondences (production  functions)  (x_.,x  ...,  x ,)  -*■ 

£ A A 01  OK 

^ ^x01 ' x02  ' •••'  xQk^  ' ^ ^ * 

ft  ft 

In  the  case  of  the  inverse  correspondence  V -*•  ]L(V)  , the 

intensity  function  Z^  is  calculated  first.  The  nodes 

A1,A2»  ...»  A^  are  taken  to  have  been  ordered  so  that  all 

intermediate  product  transfers  needed  by  the  activity  A^  are 

obtainable  from  A^,A2>  . . . , A^  ^ , i = 1,2,  ...,  k . Because 

of  this  ordering  one  can  proceed  to  node  A^_^  , imposing  as  a 

cumulative  output  demand  function  the  cumulative  inputs  required 

from  by  Ak  as  determined  by  the  computed  intensity 

function  Z^  . From  this  demand  history  one  computes  the  required 

intensity  function  for  Z ^ ^ . Then  proceed  to  node  Aj,  ^ using 

as  demand  function  the  cumulative  sum  of  the  inputs  required 

from  Ak_2  by  the  nodes  A^  ^ and  A^  as  determined  by  the 

computed  intensity  functions  Z^  , Z^_^  . In  this  fashion  the 

dynamically  interacting  intensity  functions  Zk,Z^_^,  ...,  Z^ 

can  be  calculated.  Then,  one  merely  uses  the  constraint  (2)'  as 

an  equality  to  determine  the  time  histories  xQa  , a = 1,2,  ...,  k 

of  exogenous  inputs  required  by  each  activity,  and  the  sum 
k 

x = \ x as  a time  history  of  total  inputs,  required  to 

0=1  0“  ft 

efficiently  yield  the  cumulative  output  history  V . 


21 


A A A 

In  the  calculation  of  V -*■  1L(V)  a policy  was  followed  of  producing 

as  fast  as  possible  the  outputs  required  at  each  node.  The  result- 
ing time  histories  of  intensity  levels  are  likely  to  vary  con- 
siderably from  constant  loading.  A smoothing  of  intensity  levels 
for  the  nodes  can  be  made  in  the  following  way: 


An  activity  intensity  level  Z is  critical  at  the  time  t , 

o 

if  Z ^ (t)  cannot  be  decreased  without  increasing  the  total  span 

of  time  of  the  input  history  x = £ x .By  the  policy  of 

a = l 

calculating  the  intensity  functions  Z^  (a  = 1,2,  . . . , k)  , the 

computed  functions  Z (t)  will  be  at  the  maximal  capacity  I (t) 

a a 

permitted.  If  Z is  noncritical  at  the  time  t , decrease  the 
_ a 

bound  Z a ( t ) until  Z^ (t)  becomes  critical.  By  proceeding  in 

this  way  for  all  activities  at  all  times  a new  sequence  of 

intensity  functions  Z'  is  obtained,  which  does  not  extend  the 

“ k 

total  span  of  time  of  the  input  history  x = \ x , and  yields 

a = l 

the  smoothest  loading  possible  under  such  constraint.  Further 
smoothing  is  possible  only  by  lowering  certain  intensity  capacities 
further  under  increase  of  the  total  span  of  time  of  applying 
exogenous  inputs. 


The  alterations  of  the  input  histories  obtained  by  lowering 

intensity  functions  to  criticality  are  Time  Substitutions  of  the 

calculated  required  input  histories.  The  smoothing  technique 

described  above  is  an  example  of  using  time  substitutable  input 

histories  for  economic  purposes.  One  cannot  apply  very  irregular 

loading  without  encountering  extra  cost,  and  although  the  latter 

is  not  calculated  explicitly,  the  motivation  is  clear.  All  of 

k 


the  alternative  input  histories  x = T x 

S oa 
a = l 

Z 


obtained  by  varying 


the  noncritical  intensity  functions 


without  increasing  the 

n on  the  isoquant 


total  project  span,  are  alternative  points  in  BM, 

A A 

of  the  map  set  3L(V)  , in  the  terminology  of  production^theory 
Still  other  points  on  the  isoquant  of  the  input  set  1L(V)  are 


obtained  by  decreasing  the  bounds 
extension  of  the  project  span. 


below  criticality  with 


Accumulation  of  intermediate  product  outputs  has  been  permitted, 
without  restriction.  This  may  not  be  the  case.  If  so,  a pre- 
assigned output  history  might  not  be  feasible  under  restriction 


of  inventories.  Under  such  curcumstances,  feasibility  can  be 
obtained  by  specifying  only  the  total  number  of  vessels  needed  at 
a certain  point  of  time.  Then  with  a storage-capacity-modified 
set  of  constraints  (1)'  •••  (5)'  a feasible  program  can  be  cal- 

culated and  time  substitutable  input  histories  developed  for  load 
smoothing . 

The  calculation  of  the  overall  output  correspondence  (xQ^,  • * * ' x0k^ 
P(xQ1,  ...,  Xgk)  proceeds  by  going  forward  from  activity  A^ 
to  A^  in  the  ordered  sequence  of  nodes  using  the  formulas  for 

the  intensity  functions  Z obtained  from  the  constraints  (1)  •••  (4). 

a 

Here  the  calculated  intensity  levels  Z (t)  are  either  at  their 

a 

upper  bounds  Z^  (t)  or  at  lower  levels  forced  by  the  input  histories 
x . If  an  intensity  function  Z is  not  critical  at  the  time 
t , it  may  be  lowered  to  criticality  for  load  smoothing,  without 
altering  the  output  schedule  computed  at  the  first  pass.  It  should 
be  noted  here  that  the  output  history  computed  is  merely  one 
element  of  a total  set  of  such  output  histories  obtainable  by  re- 
distribution of  the  total  input  history  x to  histories  x , 
k °a 

subject  to  £ xQa  < x . Each  such  output  history  computable  is 
a = l 

a point  in  BM™  on  the  isoquant  (production  frontier)  of  the 
output  set  2P(x)  , in  the  terminology  of  production  theory. 

A A 

Additional  points  on  the  isoquants  of  IL(V)  and  IP(x01,  . ..,  xQk) 
are  obtainable  by  altering  the  assignments  of  the  services  of 
fixed  plant  and  machinery  leading  to  the  intensity  capacity  bounds 
Z^  , so  that  alternative  capacity  bounds  are  obtained. 

The  richness  of  the  production  possibilities  in  a dynamic  model 
is  well  illustrated  here  in  this  case  of  a shipyard.  Thus  one 
can  appreciate  the  need  for  a dynamic  theory  of  production  functions 
and  the  limitations  of  our  traditional  steady  state  models. 

5.  ALGORITHMS  AND  EXAMPLE  OF  COMPUTATION 

As  an  application  of  the  foregoing  theoretical  models,  Fortran 
computer  programs  have  been  written  for  calculating  both  the 
direct  and  inverse  correspondences.  These  programs  have  been 
applied  to  activity  networks  consisting  of  up  to  22  activities 
and  time  horizons  as  large  as  150  time  units;  it  appears  that  the 
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programs  are  limited  from  handling  very  large  networks  and  time 
histories  only  by  available  computer  storage  for  variable  arrays. 
The  following  explains  the  algorithms  used,  and  the  results  of 
applying  the  programs  to  the  shipyard  example  of  Figure  1 are 
presented  as  well. 

Much  simplification  of  the  calculation  of  the  production  corre- 
spondences results  under  the  assumption  of  integral  unit  transfers 

of  intermediate  product.  That  is,  an  activity  A , with 

* ^ 

immediate  predecessor  activities  A.  , . . . , A , may  commence 

Pi  is 

i m 

production  on  output  unit  n no  earlier  than  epoch  t , where 

l t B . ) 

t = Min  (t  | l Z (a ) c 1(o)  > n , i = 1,  ...,  m>  . 

( o=l  Bi  ” ) 

In  terms  of  the  shipyard  example,  activity  A,  (outfitting)  cannot 

t h ^ 

begin  work  on  the  n1"  ship  until  activities  A2  , A-j  , A^ 

(propulsion  unit,  electric  plant,  command  control)  have  been 

t 

completed  for  the  n ship.  Such  an  assumption  is  not  seriously 
restrictive.  If  work  on  a following  activity  actually  begins  after 
a percentage  of  a given  activity  is  completed  on  an  output  unit, 
the  network  may  be  redefined,  splitting  the  given  activity  into 
two  "artificial"  activities  which  preserve  integral  output  unit 
transfers.  It  suffices  to  use  for  each  activity  a single,  scalar 
output  time  history;  for  this  purpose,  we  shall  denote  the  output 
of  activity  A^  during  [t  - l,t)  as  u^  (t)  . In  the  case 

activity  A^  has  more  than  one  immediate  follower  activity,  this 
same  time  history  u^  (t)  is  made  available  to  each  of  these 
follower  activities. 

Inverse  Correspondence 

The  algorithm  for  the  inverse  correspondence  uses  as  input  either 
the  desired  time  history  of  output  units  of  the  final  activity  or 
else  a cumulative  output  amount  at  a certain  time  horizon  (e.g. 

N ships  done  by  time  T)  to  compute  an  operations  history  for  the 


Familiar  network  terminology  is  used  throughout;  for  example, 
activity  A^  is  an  immediate  predecessor  of  Ag  means  there  is 
an  arc  connecting  A^  and  Ag  . 
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activities  of  minimum  overall  duration,  given  the  activity  intensity 
bounds  and  network  characteristics.  After  a modest  amount  of  cal- 
culation concerning  the  criticality  of  activities,  a second  itera- 
tion yields  an  operations  history  for  a set  of  minimal  activity 
intensity  peaks  (i.e.,  smoothest  activity  loading)  still  maintain- 
ing the  minimum  overall  duration  and  desired  final  output  schedule. 


By  proceeding  backwards  through  time,  the  algorithm  uses  the 
u^(t)  variables  for  accumulation  purposes  as  well.  Keeping  in 
mind  that  the  output  requirements  {u^(t)}  are  what  drive  the 
intensity  of  operation  of  activity  A^  , if  u^  (t)  exceeds  the 
output  capability  of  the  activity,  the  operation  of  A^  may  be 
set  at  capacity  and  the  unfilled  output  demand  added  onto  the 
demands  of  u^(t  - 1)  . The  algorithm  accomplishes  the  updating 

of  the  outputs  of  an  immediate  predecessor  activity  A^  by  merely 
increasing  only  outputs  { u ^ ( t ) } at  the  latest  epochs  {t}  to 
sufficiently  update  the  accumulated  output  histories. 


Given  the  order  indexing  of  the  activities  A^ , ...,  A^  , the 
algorithm  makes  use  of  the  corresponding  incidence  matrix  defining 
the  arc  network.  With  this  data,  the  algorithm  performs  the 
following  three  steps  on  the  activities  A^  , i = k,k-l,  ...,  1 . 

Step  1:  Compute  Intensity  of  Operation  for  Activity  A^ 

Given  time  horizon  T , for  t = T,T-1,  ...,  1 , set 


Zi(t)  = Min  |ui(t)/c1(t)  , Zi(t)|  . 

If  we  have  set  Z^(t)  = I^(t)  , then  set 

ui(t  - 1)  : = ui(t  - 1)  + ui(t)  - Zi(t)ci(t) 


and  set 


ui(t)  : = Z.(t)c1(t)  . 
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Step  2:  Record  Epochs  Corresponding  to  Start  of  Production  of 

Each  Output  Unit  by  Activity  A^ 

Let  N be  the  number  of  output  units  to  be  produced;  then  for 
n = 1 , . . . , N , define 


( T+1 

T. (n)  = Min  < x | £ u. (t)  > n - 1 

1 ( t=l  1 

That  is,  T^(n)  is  the  epoch  at  which  activity  A^  commences 
work  on  output  unit  n . 

Step  3:  Update  Immediate  Predecessor  Accumulated  Output  Histories 

So  As  to  Provide  Necessary  Input  to  A^ 

For  each  immediate  predecessor  activity  A ^ of  A^  (specified 
by  the  incidence  matrix) , we  require 


In  making  this 
some  such  n , 


then  define 


and  set 


test 


Ti  (n) 

£ u.(t)  > n , n = 1,  ...,  N . 

t=l  3 

for  n = 1,  N , if  one  encounters  for 

Ti(n) 

l u . (t)  < n , 
t=l  3 

T±(n) 

W = n - £ u . (t) 

t=l  3 


u^(Ti(n))  : = W + u^T^n))  . 


It  remains  to  check  if  part  or  all  of  the  required  additional 
output  (W)  is  produced  at  later  times;  if  so,  reduce  output 
at  these  times  accordingly.  Symbolically, 


for  t = T^(n)  +1,  ...,  T , set  simultaneously 
W : = Max  { 0,W  - u^ (t) } 

Uj (t)  : « Max  (0,u^(t)  - W)  . 
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Having  computed  the  intensity  of  operation  for  each  activity,  the 
schedule  of  inputs  may  then  routinely  be  calculated.  Specifically, 
for  each  activity  a , input  type  j , 

(x  (t)).  = Z (t)  aa  (t)  , t = 1,  ...,  T . 

oa  3 a ] 

Direct  Correspondence 

The  algorithm  for  the  direct  correspondence  uses  as  input  the  time 
histories  of  exogenous  inputs  xQa (t)  specified  up  to  a certain 
horizon  T , and  the  desired  accumulated  output  at  T (e.g.,  the 
number  of  ships  to  be  completed) . The  initial  iteration  produces 
an  operating  schedule  of  minimum  duration,  given  the  activity 
network  characteristics  and  intensity  bounds,  and  available  inputs. 
As  before,  a second  iteration  with  the  activity  intensity  bounds 
revised  as  to  criticality  yields  the  operations  history  with  smooth 
activity  loadings  which  still  maintain  the  minimum  overall  produc- 
tion duration. 

Given  the  order  indexing  of  activities  A^,  . . . , , the  algorithm 

makes  use  of  the  corresponding  incidence  matrix  defining  the  arc 
network.  With  this  data  the  algorithm  performs  the  following  three 
steps  on  the  activities  A^  , i = 1 , ...,  k . 

Step  1:  Determine  Epochs  at  Which  A^  is  Allowed  by  Predecessor 

Activities  to  Begin  Production  on  Each  Output  Unit 

The  algorithm  determines,  for  each  output  unit,  the  latest  epoch 
t some  immediate  predecessor  of  A^  finishes  that  output  unit; 

A^  would  then  be  allowed  to  begin  on  that  output  unit  at  epoch  t . 

For  each  immediate  predecessor  A ^ 
incidence  matrix),  for  n = 1,  ..., 

TF j (n)  = Min  j x | 

and  set 


of  A^  (specified  by  the 
N , set 

T ) 

l u . (t)  >_  n } 
t=l  3 I 


TS^(n)  = Max  { TF ^ (n)  | A^  is  an  immediate  predecessor  of  A^ } . 


For  convenience,  set 


TS^ (n)  = 0 , n = 1 , ...,N 

and  set 

TSi(N  +1)  = T + 1 , i = 1 , k . 

Step  2:  Compute  Intensity  of  Operation  for  Activity  A^  Allowed 

By  Exogenous  Inputs  and  Initial  Starting  Epoch 

For  t = 1,  ...»  TS^(l)  , set  u^ (t)  = 0 
For  t = TS^(l)  +1,  T , set 

Z . (t)  = Min  jl.(t)  , Min 

1 I 1 * 

and  set 


(i  1)  . 


(xoi(t))k 


rkl(t!  | 


ui(t)  = Zi(t)c1(t)  . 


Step  3s  Revise  Intensity  and  Ouput  of  A^  to  Reflect  Unit 
Starting  Epochs  Calculated  in  Step  1 


For  each  output  unit  n = 1,  N , the  algorithm  reduces 

outputs  of  Ai  so  that  the  cumulative  output  of  A^  does  not 
exceed  n before  epoch  TS^(n  + 1)  . Symbolically,  for 
x = TSj (n)  +1,  ...,  TSj (n  + 1)  , the  algorithm  checks  the 
cumulative  output,  and  if  for  some  such  x 


then  set 


and  set 


and  set 


l u.  (t)  >_  n , 
t=l  1 


T-l 

U.  (x)  : a n - l u.  (t) 
t=l  x 


ui(t)  = 0 , t = x + 1,  ...,  TSj  (n  + 1) 


Zf (t)  = u.(t)/c1(t)  , t = x,  ...,  TSj (n  + 1)  . 
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Example  of  Computations:  Shipbuilding  Case 

The  computer  programs,  outlined  above,  were  applied  to  data  for 
the  aggregated  network  of  shipbuilding  shown  in  Figure  1.  Although 

the  programs  in  general  will  accept  time-dependent  coefficients  

and  time-dependent  intensity  bounds  to  accommodate  learning 
effects,  in  the  following  computation  input  and  output  coefficients 
a*  (t)  and  ca(t)  were  taken  to  be  constant  through  time.  Table  1 
shows  these  coefficients  for  all  nodes  of  the  activity  network 
under  consideration. 

In  this  example,  labor  is  taken  to  be  the  limiting  input  determining 
the  operating  capacity  for  each  activity.  Thus  the  intensity 
bounds  represent  estimates  of  the  maximal  labor  input,  measured 
in  man-weeks  per  time  unit  (week) , physically  possible  for  each 
activity.  In  general,  however,  any  of  the  inputs  to  an  activity 
could  represent  the  limiting  resource  determining  an  activity's 
capacity.  Moreover,  various  activities  need  not  have  their  inten- 
sities measured  in  terms  of  the  available  input  of  the  same 
resource  or  facility  service.  Table  1 also  shows  the  capacity 
bounds  Za  (t)  , a = 1,  ...,  8 . 

Backward  Program 

With  the  above  mentioned  data,  the  inverse  (backward)  correspondence 
program  was  used  to  compute  an  operating  schedule  for  8 ships  to 
be  produced  within  a specified  planning  horizon  of  80  weeks  for 
two  final  output  schedule  requirements. 

Case  I: 


A deadline  of  80  weeks  to  complete  the  8 ships  was  specified.  An 
initial  program  run  yielded  a schedule  of  minimum  duration  beginning 
at  week  41.  The  solid  curve  of  Figure  2 shows  the  total  labor 


input  per  week  to  the  overall  network 


(1  ■-'*') 


needed  for 


project  completion.  The  uneven  loading  apparent  in  this  schedule 
necessitates  a second  calculation  to  achieve  a smoother  schedule 
without  extending  the  project  duration.  For  this  purpose,  a set 
of  minimal  constant  activity  intensity  bounds  was  computed  from 
the  first  run  and  used  in  a second  run.  Total  labor  input  per 
week  to  the  shipyard  for  the  smoothed  schedule  produced  in  the 
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second  calculation  is  shown  as  the  dashed  curve  of  Figure  2.  It 
is  to  be  noted  here  that  the  areas  under  the  solid  and  dashed  curves 
are  identical.  Figure  3 shows  the  total  material  input  (dollar 
value)  per  week  required  for  the  overall  network,  for  each  of  the 
two  runs,  while  Figure  4 shows  the  man  weeks  of  labor  per  week 
used  by  activity  A^  (ship  equipment)  for  each  of  the  two  runs. 

Case  II: 

In  this  case  the  8 ship  completion  deadlines  were  spread  through 
time,  being  demanded  8 weeks  apart  (i.e.  at  the  end  of  weeks  24,  32, 
40,  48,  56,  64,  72  and  80  respectively).  An  initial  calculation,  the 
computation  of  minimal  constant  intensity  bounds  for  smoothing,  and 
a second  calculation  using  these  bounds  were  made  as  before. 

Plotted  in  Figure  5 is  the  total  labor  input  per  week  required 
for  schedule  completion  for  each  run.  Again  the  dashed  curve  in 
this  graph  represents  the  input  flows  for  the  second  run.  Figure  6 
shows  the  total  overhead  per  week  (dollar  value)  needed  by  the 
overall  network  to  fulfill  the  above  schedule  requirement  while 
Figure  7 shows  the  man-weeks  of  labor  per  week  used  by  activity  A^ 
for  the  two  runs. 

It  is  interesting  to  note  that  in  the  second  case,  where  the  actual 
completion  times  are  demanded  from  the  shipyard  a smooth  activity 
loading  (without  extending  the  overall  project  duration)  was  only 
possible  for  A^  . In  contrast,  a complete  smooth  loading  was 
achieved  for  all  activities  in  the  first  case  where  only  the  dead- 
line for  completing  the  project  was  specified. 

Direct  (Forward)  Correspondence 

To  demonstrate  the  type  of  computation  possible  using  the  direct 
correspondence  program,  the  same  input  and  output  coefficients 
for  the  network  given  in  Table  1 were  used.  The  intensity  bounds 
Ia  , however,  were  taken  to  be  those  smoothing  bounds  developed 
from  the  first  calculation  in  each  of  the  above  two  cases  for  the 
inverse  correspondence.  Outputs  of  the  backward  program  runs  were 
also  used  as  a guideline  in  determining  a set  of  feasible  exogenous 
input  flows  for  each  activity.  Specifically,  the  maximum  input 
flow  for  each  input  type  for  each  activity  obtained  in  the  smoothed 
backward  calculation  was  used  as  a constant  input  rate  available 
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to  the  activity. 


Case  I Data: 

The  time  histories  generated  by  the  forward  program  were,  except 
for  a translation  of  starting  times,  identical  to  those  generated 
by  the  backward  program.  This  is  to  be  expected,  since  the  intensity 
bounds  and  the  available  inputs  insured  that  all  activities  had 
smooth  loading  and  no  slack  time. 

Case  II  Data: 

Although  the  overall  project  duration  and  total  inputs  used  were 
identical  to  those  generated  by  the  backward  program,  the  time 
histories  of  input  allocation  were  quite  different  from  their 
counterparts  for  the  backward  calculation.  Plotted  in  Figures  8 
and  9 are  the  total  man-weeks  per  week  used  by  the  shipyard 
necessary  for  project  completion  and  the  total  materials  used 
per  week  (in  dollar  value)  respectively. 

Only  a few  results  have  been  graphed  here  but  these  programs 
generate  complete  time  histories  of  inputs  and  outputs.  Small 
aberration  in  the  time  interconnections  of  activities  may  be  caused 
by  the  discrete  time  approach,  but  such  problems  diminish  as  the 
time  increment  is  reduced. 
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